Static triplet correlations in glass-forming liquids: A molecular dynamics study 
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We present a numerical evaluation of the three-point static correlations functions of the Kob- 
Andersen Lennard- Jones binary mixture and of its purely repulsive, Weeks-Chandler-Andersen 
variant. In the glassy regime, the two models possess a similar pair structure, yet their dynam- 
ics differ markedly. The static triplet correlation functions S*- 3 ' indicate that the local ordering is 
more pronounced in the Lennard- Jones model, an observation consistent with its slower dynamics. 
A comparison of the direct triplet correlation functions c*- 3 ' reveals that these structural differences 
are due, to a good extent, to an amplification of the small discrepancies observed at the pair level. 
We demonstrate the existence of a broad, positive peak at small wave- vectors and angles in c' 3 '. 
In this portion of fc-space, slight, systematic differences between the models are observed, revealing 
"genuine" three-body contributions to the triplet structure. The possible role of the low-fc features 
of c' 3 - 1 and the implications of our results for dynamic theories of the glass transition are discussed. 
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I. INTRODUCTION 

The origin of the vertiginous increase of structural re- 
laxation times of a viscous liquid cooled towards its glass 
transition continues to be the subject of intense research 
within the condensed matter community. In recent years, 
numerical simulations and experiments have provided ev- 
idence that the growth of relaxation times is accompa- 
nied by non-trivial dynamic correlations in the motion 
of the molecules at the nanometric scales, also known 
as dynamic heterogeneities A quantitative charac- 
terization of dynamic heterogeneities requires the evalu- 
ation of four-point dynamic correlation functions, which 
account for spatial correlations in the fluctuations of the 
dynamics [2H5| . Such a description eventually provided a 
demonstration of the role of cooperativity in supercooled 
liquids, often invoked in the past to explain the increase 
of relaxation times @ . 

Multi-point dynamic correlation functions are readily 
calculated from the particles' trajectories in computer 
simulations or in confocal microscopy measurements of 
colloidal suspensions. However, direct evaluation for 
common molecular glass-formers, such as OTP, is at 
present impossible. Notwithstanding, the four^point cor- 
relation functions can be approximated [H, 0, 0] in terms 
of three-point correlations by taking the temperature de- 
pendence of two-point dynamic functions, thereby allow- 
ing an indirect evaluation from experimental data 0, Q . 

From the theoretical side, prediction of three- 
point dynamic correlation functions has been possible 
within the so-called inhomogeneous mode-coupling the- 
ory (IMCT) [ro| . an extension of the standard mode- 
coupling theory (MCT) [ll| of the glass transition, which 
accounts for the dynamic response of a fluid to an ex- 
ternal field. The predictions of IMCT suffers from the 
same shortcomings of MCT, i.e., the dynamic correla- 
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tion volume diverges at a finite temperature; further- 
more, the predicted scaling behavior in the vicinity of 
the avoided, putative singularity is inconsistent with the 
available simulation data for four-point dynamic suscep- 
tibilities [1, [l2[- Dynamic facilitation models [l3| ex- 
plain growing dynamic correlations in terms of kinetic 
constraints on the local mobility. It can be argued, how- 
ever, that the existence of dynamic heterogeneities in 
these models is "built-in" and therefore lacks a micro- 
scopic basis. Thus, despite significant advances in the 
field, a full theoretical and experimental understanding 
of dynamical correlations beyond the bare two-point level 
is still missing. 

Along with dynamical correlations, some glass-forming 
liquids appear to develop non-trivial static correlations 
upon cooling below the onset temperature. Indeed, sev- 
eral recent studies have successfully identified static cor- 
relation length scales and studied their temperature de- 
pendence [lj-[lj|- The rationale behind these studies is 
that not only the dynamics, but also the local structure 
is heterogeneous. Such static heterogeneities may corre- 
spond to domains formed by locally preferred structures, 
as those envisaged in frustration-based theories of the 
glass transition [l9l - [2l| . or to the appearance of some 
more exotic amorphous order [Til [23 - 124^ . 

Simulation studies of both simple and realistic mod- 
els of metallic glass-formers have given concrete shapes 
to locally preferred structures j25l - l35| . e.g., icosahedra or 
more general polyhedral structures. Almost invariably, 
such static heterogeneities are found to be connected to 
the dynamic properties, such as dynamic heterogeneities 
or frag ilit y . Simulation [3(1 and experimental stud- 
ies [3|| |39[ of bidispcrsc or polydisperse soft particles in 
two dimensions reveal how the competition between hex- 
atic order and defective local structures impact on the 
dynamics. Further evidence of the connection between 
static and dynamic length scales [4(| is provided by sim- 
ulations of Lennard-Jones particles confined on the hy- 
perbolic plane. On the other hand, a related, a recent 
study [16| points towards a marginal role of static corre- 
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lations in hard sphere models, indicating that the exis- 
tence of a link between structure and dynamics may be 
system-specific. 

Properties used to characterize static heterogeneities, 
such as bond-order parameters, are often tailored to de- 
tect specific types of local order. Furthermore, differ- 
ent geometrical constructions, e.g., Voronoi tesselations 
and common neighbor analysis, may lead in some cases 
to inconsistent results |4l| ■ Therefore, a more general 
and "order-agnostic" framework to evaluate static corre- 
lations is desirable. A recent line of research builds on the 
assumption that the relevant static correlations in glassy 
systems are of high order, i.e., well beyond the pair level. 
Indeed, Voronoi and common neighbor analysis typically 
involve correlations between several particles at a time 
within the first neighbor shell. One general method to 
account for high order static correlations in a liquid is 
to evaluate the so-called "point-to-set" correlations 
This technique requires pinning a fraction of the particles 
and evaluating the influence of the external field of the 
immobile particles on the structure of the liquid [l8| . At 
least for one model glass-former, the popular Lennard- 
Jones Kob- Andersen (KA) mixture, point-to-set correla- 
tions have been shown to possess a non-trivial connection 
to the dynamics [l?}- Specifically, Hocky ct al. have ex- 
plained the differences in dynamic behavior between two 
variants of the KA model (with and without attractions) 
in terms of different rates of growth of the point-to-set 
lengths. Their results are consistent with a previous anal- 
ysis of the temperature variations of the local structure 
of these models [42j . 

Here we take a different point of view. Instead of 
calculating many-body correlations functions, without 
knowing exactly which order is involved, we ask ourselves 
whether it is possible to extract relevant information on 
the glassy dynamics already from few-body correlation 
functions. In our previous work [42J it was found that 
the distribution of angles formed by nearest neighbors is 
more sensitive than the pair correlation functions as an 
indicator of the differences between attractive and repul- 
sive variants of the KA model. In the present work, we 
perform a systematic evaluation of the three-point static 
correlations of these two models. We find that a sub- 
stantial part of the difference in the triplet correlations 
arise from an amplification of the small differences in the 
pair correlations within the first coordination shell. The 
direct triplet correlation functions of the models differ 
slightly, but systematically, at larger length scales. Fi- 
nally, we briefly discuss the different role of direct triplet 
correlations in close-packed and tetrahedral network liq- 
uids. 

The manuscript is organized as follows. In Section |TT] 
we briefly describe the models and present our results. 
In Section Mil we discuss the implications for dynamic 
theories of the glass transition and give our conclusions. 



II. MODELS AND RESULTS 

We report molecular dynamics (MD) simulations for 
two versions of the Kob-Andcrscn (KA) binary mix- 
ture [43|: in the original model, particles interact via 
the Lennard- Jones (LJ) potential 
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where a, j3 = 1, 2 is a species index. The values of the pa- 
rameters can be found in the original paper |43l . The po- 
tentials are cut and quadratically shifted [44| at 2.5a a p. 
We also study the so-called Weeks-Chandler-Andersen 
(WCA) variant of the model 45J, in which the poten- 
tials u a p(r) in Eq. (TTJ) are truncated and shifted at the 
minimum and are there fore purely repulsive. The pa- 
rameters (Tii, en, and ^/miCjj/en constitute our units 
of distance, energy, and time, respectively. We simulate 
systems of N = 1000 particles in a cubic box with peri- 
odic boundary conditions at a constant number density 
p = N/V = 1.2. All simulations were carried out in the 
NVT ensemble using the Nose-Poincare thermostat [4(| 
with a mass parameter (5=5.0. Some larger samples of 
N = 8000 have also been simulated to check the behavior 
of static correlations at large length scales. 

It has been pointed out that LJ and WCA models ap- 
pear similar at the level of static pair correlations 14 711 
but display subtle differences in the local structure [42j |. 
Around its putative mode-coupling critical temperature 
(Imct ~ 0.435), the LJ model is around two orders 
of magnitude slower in relaxation times than the WCA 
model. This stark difference in dynamic behavior has 
been related to the different thermal rates of growth 
of clusters formed by locally preferred structures [isT ]. 
as well as to the different growth rates of the point-to- 
set lengths [I?}- Due to the strong similarity of their 
two-point structure, LJ and WCA models offer an ideal 
benchmark for testing the role of static multi-point cor- 
relations beyond the pair level. In the following, we will 
evaluate the three-point static correlation functions, with 
the additional aim of providing guidelines how to im- 
prove microscopic theories of the glass transition based 
on few-point static correlations, such as MCT or replica 
theory |49j . 

In multi-component systems, extracting a sensible 
physical picture from the triplet correlations is difficult 
due to the large number of parameters' combinations to 
consider [Hoi l5l| . Binary mixtures are characterized by 
8 independent triplet functions, each one being a func- 
tion of, at least, three scalar variables. In order to tackle 
the complexity of parameters space and to simplify the 
presentation, we have chosen to focus only on selected 
triplets, a/?7, of species. Our choice is motivated as fol- 
lows. We will report results for 111 correlations, because 
the large particles (a = 1) constitute the majority species 
in the KA model; the corresponding correlations are thus 
expected to impact substantially the outcome of dynamic 
theories, e.g., through the memory kernel in MCT. In 
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FIG. 1. Triplet correlation functions (a) gu{(r, r, cos (9)) and 
(b) gif^r, r, cos (6»)) for the LJ model (full lines) and the WCA 
model (dashed line) at T = 0.45. 



addition, we study 121 correlations, because structural 
motifs in binary alloys are known to be centered around 
the minority (solute) component [52ll54j . i.e., the small 
particles in the KA model; these correlations will thus be 
more sensitive to the details of the local structure. We 
have carefully inspected all species triplets to check the 
robustness of our conclusions. 

The physical content of the triplet correlations is most 

easily grasped in real space. Let us therefore start our 

presentation with a brief analysis of the triplet Corre- 
al -* 

lations g^p^r, r), where r and r are the vectors con- 
necting a central particle of species /3 with two parti- 
cles of species a and 7 at distances r and r*, respec- 

tively. In an isotropic, homogeneous liquid, g a/ 3 7 (r, r) = 
(3) 

5a/3 7 ( r ' r '> cos (0))> wnere is the angle formed by the two 
vectors r and f. In Fig. [T] we show g-^ and g^i for the 
two KA models at T = 0.45, i.e., close to the estimated 
Imct of the LJ model. For simplicity, we only consider 
here distances along the diagonal, r — r' . The differences 
between the LJ and WCA models appear more evident in 
the 121 than in the 111 triplet correlations. The ampli- 

(3) 

tude of the peaks of g\%\ are in fact systematically larger 
in the LJ model than in the WCA. The more pronounced 
local ordering of the LJ model at low temperature is 
consistent with its slower dynamics. By contrast, the 



two models display rather similar correlations between 

triplets of large particles. In particular, there is barely 

(3) 

any difference in the amplitude of the main peak of g^n, 
which arises from angles « 60 between three neighboring 
large particles. The peak is only shifted to slightly larger 
distances in the LJ model. This demonstrates that the 
structural differences between LJ and WCA models are 
mostly due to the subtle ordering of the large particles 
around the small ones. Thus, static properties involving 
correlations between unlike species are best suited to de- 
tect structural changes in the KA model, which are due 
to the appearance of bicapped prismatic structures [42[ 
centered around small particles. As recent work indi- 
cates !5j|, these observations may prove useful in the 
analysis of the low-activity, glassy configurations gener- 
ated by biasing the particle mobility [5(1 [57| , or of the 
ultra-stable glasses obtained by vapor deposition tech- 
niques [58j . Finally, we point out that the bond-angle 
correlation function Di2i(9) presented in Ref. [42j can 

be obtained as the integral of <?i 21 (f\ r) over the first 
shell of neighbors. 

From the point of view of microscopic dynamic the- 
ories, correlations are most conveniently represented in 
Fourier space. In the following, we will therefore focus 
on the Fourier transform of the triplet correlation func- 
tion 



^7(^9) = jr(p a (q)P0(k)p*(k + q)} 



1 



iV 



(2) 



where a,/3, and 7 are indices of species and p a (k) is the 
Fourier transform of the microscopic density p a {k) = 
Sj^i exp(z/c • fj). Wc will often refer to the ensemble 
of all possible triplet correlations simply as S^ 3 \ Under 
hypothesis of homogeneity and isotropy, S^ 3 " 1 is a func- 
tion of (k, g, cos#), where 9 is the angle formed by k and 
q. For the sake of simplicity wc will mostly concentrate 
ourselves on slices in Fourier space corresponding to diag- 
onal wave- vectors, k = q. Our conclusions are confirmed 
by an analysis of the triplet correlations as a function of 
k and q at fixed angle 6. For completeness, a brief discus- 
sion of triplet correlations for off-diagonal wave-vectors 
will be included at the end of the section. 

Figure [2] shows the triplet correlation functions 
s[ 3 1 > 1 (k, k, cos(0)) and S^ 1 (k,k,cos(6)) evaluated at T = 
0.45. Triplet correlations between large particles display 
two prominent peaks along the cut k = q = 7, for angles 
cos(0) w —0.5 and cos(6') w 0.45, respectively (see also 
Fig-El)- The dominant wave- vector (k* = 7) is also the lo- 
cation of the first peak in the structure factor Sn(k). To 
pinpoint the physical origin of these "preferred" angles, 
we tentatively compare the peaks' intensities to those in 
real space (Fig. []}. For instance, the largest peak at 

(3) 

wide angles in could be attributed to the preferred 
angle of 9 « 60 between three neighboring large parti- 

(3) 

cles. In the case of <Si 21 , the negative minimum around 
k = 6 arises in correspondence to the first minimum of 
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FIG. 2. Triplet correlation functions (a) s[f\ (k, k, cos (#)) 
and (b) S^iCfc, k, cos (#)) for the LJ model (full lines) and 
the WCA model (dashed line) at T = 0.45. 



the structure factor Si2(k). Inspection of as a func- 

(3) 

tion of q and k at fixed angle 0, shows that S^H is sym- 
metric in q and fc, as expected, and takes its absolute 
maximum on the diagonal. The absolute maximum of 
is attained for slightly off-diagonal wave- vectors. 

As expected, the systematic differences visible in S^i 
confirm the more pronounced ordering of the LJ model. 
In particular, the amplitude (in absolute values) of the 
peaks of S^i is systematically larger in the LJ than in 

(3) 

the WCA model. Differences in S^H are only visible at 
small angles; by contrast, barely any difference is ob- 
served around the main peak. These features are further 
highlighted in Fig. [3J which shows the angular dependen- 
cies of the s[f[ and S^i a * fixed wave- vector q = k = k* . 
The different choices of k* for and S^i reflect the 
locations of the meain peaks of Sn and S12, respectively. 

We found that triplet correlations are fairly sensitive 
to temperature variations — much more than pair corre- 
lations. The locations of the peaks are only marginally 
affected by temperature variations, but their amplitudes 
(in absolute value) show a substantial increase upon cool- 
ing. To quantify the latter effect, we identify the abso- 
lute maximum of over all wave-vectors k and q and 
track the evolution of its height as a function of 

T, see Fig. @] We inspected the location of the absolute 
maxima to make sure the peaks evolve smoothly, i.e., the 
geometrical meaning of the peak is consistent across tern- 
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model (full lines) and the WCA model (dashed line) at T — 
0.45. 



peratures. For comparison, we also include the heights 
Sn and S12 of the main peaks of the structure factors 
Sii(fc) and Si2(k). The temperature dependence of 
is stronger than that of S. The variation of the absolute 
maxima of S^ 3 ' covers a factor 6-7 (depending on species 
triplet) going from high (T = 4) to low (T = 0.45) tem- 
perature. By contrast, the peaks of the structure factors 
increase by a more modest factor two. Further, it can be 
seen that, below the onset temperature (T w 1.0) of the 
LJ model, S\ 21 increases at a higher thermal rate in the 
LJ model than in the WCA. As it may be expected, how- 
ever, the difference is not as marked as the one revealed 
by the analysis of the locally preferred structures [42| or 
by higher order correlations [17] . Of course, our simpli- 
fied analysis carries no information on how the spatial ex- 
tension of triplet correlations changes with temperature. 
Whether a growing static length can be extracted from 
three-point correlations is an issue that deserves further 
investigation. 

The above results demonstrate that static correlations 
become more sensitive to temperature variations with in- 
creasing order. This, in turn, makes it easier to disentan- 
gle the behaviors of the LJ and WCA models. Careful 
inspection of the static structure factors (or the radial 
distribution functions) reveals, however, that small, sys- 
tematic differences can be discerned already at the two- 
body level. It is thus natural to ask whether the differ- 
ences observed in the three-point correlations arise from 
an "amplification" of those already visible at the two- 
body level, or reflect some "genuine" many-body effects. 
A similar issue has been recently discussed in connection 
to point-to-set correlations (59j . 
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At the level of three-point static correlations, the above 
question can be phrased precisely in the language of liq- 
uid state theory: to what extent the convolution approx- 
imation accounts for the differences in S 1 ' 3 ' between the 
two models? In fact, the convolution approximation fac- 
torizes the triplet correlations into products of structure 
factors, i.e., two-point functions; for a multi-component 
system it reads 



Np 2 [S a a'(k)S m (q)S jY (\k + q\) 

a' /9'7' 

$a<i3'Sa'~/'S~/'l3<] (3) 



Note that, in contrast to the simple one-component case, 
for which 



'(k,q) = S(k)S(q)S(k + q), 



(3) 

the triplet correlation, say, S^H receives now contribu- 
tions from Sip, with arbitrary /3. The direct triplet cor- 
relation functions for a multi-component system are then 
denned by the linear system [6(| 



(4) 



s { X(k,q) = N P < 



„2 \ " 



S aa '{k)Spp-(q)S 11 >{\k + q\) 



[5 a >i3'5 a 'p'5a'P</Xo 



P^&q)) (5) 



This somewhat complex expression hides the physical 
meaning of the direct triplet correlation, which is evident 
in the one-component case: 

S^(k,q) = S(k)S(q)S(k + q)[l + p 2 c^(k,q)] . (6) 

c (3) j s proportional to the relative deviation of the actual 
triplet correlation from the convolution approximation 
5 conv . Indeed, the convolution approximations, Eq. [3] 




FIG. 5. Triplet correlation functions S {3) (k, k, cos (0)) (full 
lines) and the convolution approximation, Eq.[3] (dotted line) 
for the LJ model at T = 0.45. (a) 111 correlations; (b) 121 
correlations. 



and Eq. 0] arc obtained by setting = in Eq. (0) and 
([6]) , respectively. It is in this sense that c 1 - 3 ' is considered 
as a measure of "genuine" three-point contributions to 
the triplet correlations. 

To test of the validity of the convolution approxima- 
tion, we first compare S^ 3 ) and S conv . From Fig. [5] we 
see that the convolution approximation works fairly well 
along the diagonal k = q. In particular, the major peak 

(3) 

at wide angles in SlU is accurately reproduced, and so is 
its temperature variation. Differences arc visible at small 
angles, where the profiles at constant k obtained by con- 
volution appear flat and fail to reproduce the oscillations 
in the angular dependence. Nonetheless, an inspection of 
the triplet correlations for off-diagonal vectors, as well as 
for the WCA model, confirms that that the convolution 
approximation works overall reasonably well. This result 
may be expected in dense, hot liquids, but is non-trivial 
close to the estimated MCT critical temperature. These 
observations can be made more quantitative by analyzing 
the absolute errors introduced the convolution approxi- 
mation (not shown here) . The absolute errors are indeed 
largest at k w k* and small angles, and become negligi- 
ble (typically < 10~ 2 ) at wave- vectors smaller than k*. 
Still the relative error at small k can be very large, as 
discussed in the following. 

To what extent the differences in between the 
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two models arise from correlations not visible at the pair 
level? To answer this question we evaluate the direct 
triplet correlations. Let us first inspect of the behavior 




of and c^i at wave- vectors q = k ps k* (sec Fig. [BJ. fig. 7. Direct triplet correlation functions (a) 



,(3) 



At these wave-vectors, c 1 - 3 ' resembles the familiar shape 
known from earlier studies on dense liquids with hard 
core interactions [6l[. Strikingly, c^\(fc*, k*, cos (6)) 
matches exactly in the two models; minor differences in 
c 121 (k* , k* , cos (9)) are visible at small angles, where the 
convolution approximation slightly underestimates cor- 
relations in the LJ model. A direct comparison of S conv 
in the two models confirm these observations. We con- 
clude that, for diagonal wave- vectors around k* , the con- 
volution approximation introduces systematic errors in 
the triplet structure, but these errors are similar in the 
two models. A possible physical interpretation is the fol- 
lowing: c^(k*, k*, cos (#)) is sensitive to the geometri- 
cal details of local ordering in the first neighbors shell. 
Since LJ and WCA models possess the same type of lo- 
cally preferred structure [42 1, the direct triplet correla- 
tions around k* look very similar. The difference lies in 
the strength of the correlations, which is already well ac- 
counted for by the convolution approximation, and their 
spatial extension, which might appear at smaller wave- 
vectors. 

Inspection of C3 at wave-vectors lower than k* reveals 
indeed small, but systematic discrepancies between the 
models (see Fig.[7]). Both ctf^k, k, 9) and c{2i(k, k, 9) de- 
velop, at small k and small angles, a large positive peak, 
which is slightly higher in the LJ than in the WCA model. 
We found that the amplitude of this peak increases by 
decreasing temperature. Further, an investigation of the 



Cjjj (k, k, cos (0)) and (b) c^ 1 (k, k, cos (0)) for the LJ 
model (full lines) and the WCA model (dashed line) at 
T — 0.45. The range on the z-axis has been limited from 
below for clarity. 



direct correlations for off-diagonal wave- vectors confirms 
that discrepancies in c^ 3 \ whenever present below k* , im- 
ply a more pronounced local ordering in LJ model than 
in the WCA one; see Fig. [5] for representative results at 
selected small angles. To the best of our knowledge, the 
existence of such a broad maximum in the direct triplet 
correlation of a fluid has never been explicitly demon- 
strated in the literature. A positive peak in c^ 3 ' at small 
k and small angles was reported by Jorge et al. [5 1 for 
an hard-sphere mixture. The amplitude of this peak was, 
however, very small, probably due to the relatively low 
packing fraction of the model employed in that study. 
We performed additional simulations for the hard sphere 
model of Ref. [5l| and confirmed the presence of this 
small maximum. At even smaller k, the direct correla- 
tion functions eventually become very large and negative 
for all species triplets, in agreement with earlier studies 
on hard spheres [5l|, |6l|. In this portion of fc-space, a 
clear-cut comparison between the models is made diffi- 
cult by the large statistical uncertainties. 

The above results indicate that at intermediate length 
scales there arc small contributions to the triplet corrcla- 
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9 = 60 6 = 46 6= 18 

FIG. 8. Direct triplet correlation functions c^-^fe,^ cos (#)) and for the LJ model (full lines) and the WCA model (dashed 
line) at selected angles (as indicated in the labels) at T = 0.45. 



tions, ignored by the convolution approximation, which 
reinforce the structural difference between the models. 
Whether the low-fc structure of has any physical rel- 
evance is unclear at present. It can be seen from Fig. [2] 
that S^ 3 ) is negligible in the corresponding portion of k 
and 8, and so is the absolute deviation from S conv . In 
its standard formulation, MCT gives a strong weight to 
wave-vectors close to the main peaks of the static struc- 
ture factors, due to the quadratic nature of the mem- 
ory kernel. Thus, one may expect contributions at small 
wave-vectors to be suppressed, effectively washing out 
any difference between LJ and WCA models beyond the 
pair level. In close-packed liquids, inclusion of affects 
marginally the solution of the MCT equations, while a 
stronger effect is found in network liquids (62j . This has 
been tentatively attributed in Ref. [62| to positive con- 
tributions to the memory kernel arising from off-diagonal 
terms in c^ 3 \ To address this issue, we also studied the 
triplet static correlations of a simple model of tetrahc- 
dral network liquid [(33[ . Preliminary results indicate the 
existence of positive peaks at small k, similar to those 
observed for LJ and WCA models, but not for all species 
triplets. In the case of the network model, the small-fc 
peaks of cS 3 ^ appear in the range of wave-vectors corre- 
sponding to the first sharp diffraction peak of the struc- 
ture factors. This portion of fc-space is thus still sig- 
nificantly weighted by the MCT memory kernel, which 
might explain the impact of c' 3 ^ on the solution of the 
MCT equations for network liquid. Work to test this 
hypothesis is in progress. 



III. CONCLUSIONS 

The comparison between LJ and WCA models has 
been taken by Berthier and Tarjus [13, I64l466j as a 
paradigmatic exemple of the limitations of microscopic 
theories of the glass transition based on the sole pair 
structure. High-order static correlations have therefore 
been invoked as a means to disentangle the dynamic be- 
haviors of L J and WCA models [l?], EJ . In this work, we 
have analyzed the three-point static correlations of these 
models. We found that the triplet correlations S^ 3 \ even 
when approximated by convolution, are much more sen- 
sitive to temperature variations than two-point correla- 
tions. Further, we showed that differences in the triplet 
structure of the LJ and WCA models arise, to a good 
extent, from an amplification of the small discrepancies 
visible at the pair level. The direct triplet correlation 
functions are in fact surprisingly similar for wave- 
vectors close to main peak of the structure factor. Slight, 
systematic differences in are observed at smaller k, 
where a broad, positive peak develops for certain species 
triplets. Whether the low-fc structure of may expose 
a static length scale relevant for the glass transition is a 
question deserving further study. 

The results of this study help clarifying the role of the 
approximations involved in dynamic theories of the glass 
transition based on few-point static correlations, such as 
MCT. One should bare in mind, for instance, that both 
three- and four-point static correlations are indeed taken 



8 



into account by MCT, although in a factorizcd form. Us- 
ing the pair structure factors as input, a solution of the 
MCT equation (6f| only partly accounts for the dynamic 
differences between the LJ and WCA models. Given the 
strong similarity of their and the negligible effect of 
c^ 3 ) on the MCT solutions for close-packed liquids j6^ |. 
the deficiencies of the theory will not be cured by inclu- 
sion of the direct triplet correlation functions. Alterna- 
tive approximation schemes for the memory kernel [67j . 
or different microscopic approaches, such as replica the- 
ory [G^, may be able to amplify the small- k features 
of and thus prove more effective than the standard 
MCT in disentangling the dynamics of the two models. 
On a more general note, our work indicates that system- 



atic inclusion of factorized static and dynamic correla- 
tions beyond the pair level, such as in generalized ver- 
sions of MCT [o7l l69| , is a promising line of research to 
improve constructively existing dynamic theories of the 
glass transition. 
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